Temperature controlled Levy nights of minority carriers in photoexcited bulk ?i-InP 
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We study the spatial distribution of minority carriers arising from their anomalous photon-assisted 
diffusion upon photo-excitation at an edge of n-lnP slab for temperatures ranging from 300K to 
78K. The experiment provides a realization of the "Levy flight" random walk of holes, in which the 
Levy distribution index 7 is controlled by the temperature. We show that the variation 7(T) is close 
to that predicted earlier on the basis of the assumed quasi-equilibrium (van Roosbroek-Shockley) 
intrinsic emission spectrum, 7=1 — A/kT, where A(T) is the Urbach tailing parameter of the 
absorption spectra. The decreasing 7 at lower temperatures results in a giant enhancement in the 
spread of holes - over distances exceeding 1 cm from the region of photo-excitation. 

PACS numbers: 78.30.Fs,78.55.Cr,78.60.Lc,74.62.En 



Random walk of the Levy flight kind arises in many 
diverse fields, including hydrodynamics biology 
financial markets Q, earth science etc. It is very im- 
portant whenever the photon transport is mediated by 
multiple photon absorption-reemission processes in gases 
in particular in astrophysics Q. Recently, the Levy 
flight was discovered in photon-assisted transport of mi- 
nority carriers in semiconductors Q. 

In the normal transport process, the stationary distri- 
bution of minority carriers produced by optical excitation 
in semiconductors decays exponentially from the excita- 
tion area and is characterized by a micron-scale diffusion 
length I. However, in moderately-doped direct-gap semi- 
conductors with high radiative efficiency, the minority- 
carrier transport is strongly modified by the "photon re- 
cycling" process (repeated radiative recombination and 
reabsorption of emerging photons). The distribution of 
steps in this photon-mediated random walk of minority 
carriers is defined by the photon reabsorption probability 
0] and the single-step probability distribution V(x) can 
be calculated from the interband absorption and radia- 
tive recombination spectra Q. 

At the red edge of the spectra in the Urbach-tail region 
the reabsorption length becomes anomalously large. As 
a result, V(x) asymptotically obeys the power law 0] 



V(x) ~ 1/x 1 



+1 



(1) 



with < 7 < 1. Since the distribution ([1} is "heavy- 
tailed," the diffusivity [conventionally denned through 
the second moment of V(x)] diverges and the random 
walk is governed by rare but large steps. This kind of 
transport is known as the Levy flight Tol- 12|. 

The resultant stationary spatial distribution is lim- 
ited by: (i) loss of photons in the free-carrier absorption 
(FCA) and (ii) loss of minority carriers in nonradiative 
recombination. In n-doped InP the room-temperature 
FCA coefficient is a F cA [cm" 1 ] « 1.3 x 10~ 18 iVn_[cm- 3 ] 
and decreases further at lower temperatures [13j, 
Thus, for moderately doped samples this mechanism 
of loss is ineffective. The relative rates of recombina- 



tion are characterized by the radiative efficiency 77 = 
T nr /(T n r + T rad ) (where t"* and t^ 1 are, respectively, 
the radiative and nonradiative recombination rates) . The 
average number of steps is given by the recycling factor 
$ = 77/(1 — 77), which is typically $ 3> 1. 

In the previous study Q , we directly observed a power- 
law decay of the hole concentration, characteristic of the 
Levy flight, function of the distance from the ex- 
citation — up to several millimeters — for differently 
doped samples at 300 K. Here we demonstrate that the 
index 7 of the Levy distribution can be controlled by 
varying temperature. The decrease of 7 at lower T pro- 
duces further enhancement in the hole spread, reaching 
centimeter-scale distances. 

The temperature dependent luminescence was studied 
in the same geometry as earlier Q , using n-InP slab [lH 
of size 20 mm by 7 mm and thickness d=350 /^m, cor- 
responding to x, y, and z directions, respectively. The 
moderate doping level Ad = 3 x 10 17 cm~ 3 gives op- 
timum radiative efficiency at room temperature A 
808-nm laser beam was focused on the 7-mm edge by 
a cylindrical lens producing a uniform excitation along 
the edge. The excitation energy E = 1.53 eV was well 
above InP bandgap (E g « 1.35 eV) and generated holes 
only in a submicron layer near the edge. We registered 
the intensity and spectra of luminescence emitted from 
the broadside as a function of the distance x from the 
photo-excited edge. 

The luminescence spectra retain identical shape in the 
whole range x. This important observation remains valid 
for all temperatures studied. The edge-excited photolu- 
minescence spectrum for T = 296 K, observed at x=l 
mm, is shown in Fig. [T] together with the calculated 
spectrum. Also shown is the luminescence spectrum, ob- 
served from the broadside in reflection geometry ja, Q ■ 
We note that, compared to the reflection spectrum, the 
edge-excitation spectrum is noticeably shifted to the red 
side, indicating substantial filtering [7[. 

The absorption spectrum a(E) of the sample is shown 
in Fig. [T]by the dash-dotted line. In the whole range of 
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FIG. 1. (Color online) Luminescence spectra of n-doped 
InP sample. Photo-luminescence spectra (1) for edge exci- 
tation and front-side observation; experimental data (circles) 
are practically indistinguishable from the calculation (dashed 
line) that allows for spectral filtering. Reflection luminescence 
spectra from the front side is shown by the solid line (2). The 
van Roosbroek-Shockley intrinsic emission spectrum, Eq.© 
is shown by the dashed line (3). The dash-dotted line (4) 
corresponds to the absorption coefficient. 



temperatures studied, the red wing of a(E) exhibits an 
exponential Urbach behavior (l3j . 



a.i{E) = ao exp 



E - E„ 



A(T) 



E < E a 



(2) 



At E > E g , the Ur- 



extending down to a(E) ~ apcA- 
bach exponent saturates. Note that that the lumines- 
cence spectra observed with the edge excitation are fully 
in the Urbach tail region of the absorption spectrum - 
in contrast with the intrinsic emission spectrum. 

Assuming a quasi-equilibrium energy distribution of 
holes, the intrinsic emission spectrum is given by the well- 
known van Roosbroek-Shockley (VRS) relation [l6| . 



S VRS (E)~ai(E)E 3 exp 



E 
kT 



(3) 



The VRS spectrum is also shown in Fig. [TJ It has a 
maximum above E g . Owing to the rapid energy relax- 
ation of holes, the intrinsic emission spectrum in n-InP 
is well described by Eq. (O, both at room and lower 
temperatures, with a noticeable deviation in shape only 
at T < 78 K. The observed red-shifted spectrum S(E) 
is shaped by reabsorption on its way out of the sample, 
S(E) = Svrs x F(E), where F(E) is the spectral filter- 
ing function F(E) = F X {E) T(E) that depends on the 
hole distribution p(z) across the wafer and is affected by 
reflections from the sample surfaces. The one-pass filter- 



ing function F\ {E) is given by 0, Q 

F X {E)= [ p{z)cxp[-a(E)z]dz 
Jo 



(4) 



The factor T(E) = [1 - i?cxp(-a(S)d)]- 1 accounts for 
the multiple surface reflections; R sa 0.33. Note that due 
to the high index contrast the radiation escape cone is 
narrow, i.e. the outgoing radiation propagates close to 
the normal direction to the surface. Effects of multiple 
reflections are noticeable only in the far red wing of the 
spectrum. For the edge excitation, the hole distribution 
across the sample is nearly homogeneous. Therefore, the 
typical radiation escape distance in this geometry is much 
larger and produces considerably larger red shift of the 
spectrum than for the reflection geometry, where holes 
are generated near the surface. The filtered spectra, cal- 
culated using the VRS relation, are shown in Fig. [T]by 
the dashed line (1), demonstrating excellent agreement 
with the experiment 17 1. 

The luminescence intensity distributions I(x) were ob- 
tained by scanning the CCD image along a line parallel 
to x, as in the earlier experiment Q- The proportion- 
ality p(x) cx I{x) is supported by the good agreement 
between the calculated and observed luminescence spec- 
tra for edge excitation. 

The resulting distributions for temperatures ranging 
from 300 K to 78 K are shown in Fig. [2l One can see a 
huge enhancement of the hole spread that extends over 
1 cm at temperatures below 200 K. 

The power law p{x) ~ 1 /x l+ ~' is clearly observed for all 
temperatures at distances x > 0.5 mm. This is best seen 
on the log-log scale in the inset of Fig. [2j The power 
law is in clear contrast to an exponential decay p(x) ~ 
exp(— x/l) expected for a normal diffusion of holes, even 
accounting for any photon-assisted enhancement of the 
diffusion length I. 

The observed hole distribution p(x) for all temper- 
atures follows the calculated distribution for the Levy 
flight random walk governed by the single-step probabil- 
ity |T]). Details of this calculation have been extensively 



discussed @, H, 18,[19|. The V(x) itself is evaluated as 
the photon reabsorption probability @, S] ; it is fully de- 
termined by the absorption and the emission spectra. It 
can be calculated numerically in the full range of dis- 
tances. At large distances (x 3> a$ 1 , 1) a tangible contri- 
bution to V[x) results only from the exponentially vary- 
ing Urbach part of the absorption spectra. This region 
also correspond to the exponentially decaying red wing 
of the intrinsic emission spectrum. Assuming the intrin- 
sic spectrum in the VRS form, we find [1, Q an analytic 
expression for index 7 of the single-step distribution 



1VRS 



kT 



(5) 



An analytic expression for the hole distribution p(x) 



3 




2000 4000 6000 8000 

Distance, u.m 



10000 12000 



FIG. 2. (Color online) Distribution of the luminescence inten- 
sity I(x) and the hole concentration p(x) oc I(x) for different 
temperatures T. Dashed lines correspond to the power-law 
([T]) fitting with an index 7. The inset shows the same distri- 
butions in log-log scale. 
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FIG. 3. (Color online) Temperature dependence of the pa- 
rameters of Levy distribution. Index 7 shown by squares is 
estimated from the long-distance asymptotics of the measured 
p(x); also shown (triangles) are the theoretical values (0. The 
excitation front distance Xf is estimated from maximum cur- 
vature in the log-log plots of Fig. [2] Lines are guides to the 
eye. 




can be derived ]2 
proximation [(| [2 

p(x) 



in the so-called "longest flight" ap- 
and is of the form 



x l+7[l + (xf/x) 



712 



(6) 



where c is a normalization constant. Distribution (|6]) 
provides a good approximation to the exact solution in 
the entire range of x. Asymptotically (at x xf) it 
reproduces the single-step probability ([T]), exhibiting a 
power-law decay with index 7. At short distances, x <C 
Xf, it gives a weaker decay, p(x) cx x 7_1 , in agreement 
with the exact solution 0, Efl ■ The distance x / gives an 
estimate of the spread of the excitation "front" of p(x) 
— beyond which the holes appear predominantly in one 
step. The transition between the short- and long-distance 
asymptotics is clearly seen for all T in the inset to Fig. 
with x f being the point of maximum curvature in the 
log-log plots. 

Wc have estimated 7 through the slope of log-log varia- 
tion at far distances and also x / via the position of maxi- 
mum curvature oip(x). The temperature variation of ob- 
tained values of 7 and Xf is shown in Fig. [3] Also shown 
are theoretical values of the index jvrs calculated §5§ 
by assuming the intrinsic emission spectrum in the VRS 
form ((3]). We see that the estimation of 7 through the 
Urbach tailing parameter gives a good agreement with 
experiment. 

Position of the excitation front can also be estimated 
through the recycling factor, viz. Xf = xo<& ' 7 , 
where xq « 0.2 (im, is the length scaling parameter of 
V(x). The observed huge increase of the excitation front 
distances, 1/ > u P to the mm-range at low temper- 



atures is another hallmark of Levy flight, which cannot 
be attributed to any other mechanism of hole transport. 
We note however, that for $ « 100 and observed j(T) 
one should expect even higher increase of Xf. Therefore, 
our experiments suggest a decreasing recycling factor at 
lower T . Apparently, the well-known increase of the ra- 
diative recombination rate at low T is accompanied by an 
even faster growth rate of non-radiative recombination, 
as found for several recombination models with vanishing 
activation barrier 22h24|. 

From a practical perspective, the anomalous transport 
of minority carriers should be important in all semicon- 
ductor devices with high radiative efficiency Q. How- 
ever, details of the transport and the temperature effects 
depend on the absorption spectrum in the tailing region 
and are not universal, e.g. often deviate from the Urbach 
shape. Thus, some semiconductors feature Gaussian tails 
[2^ . We have calculated the index 7 for Gaussian tails 
and found it to increase with decreasing temperature, re- 
maining within a narrow range (0.85 < 7 < 1.05) close 
to unity for a wide temperature variation. Such semicon- 
ductors may serve as an experimental model for study- 
ing the transition from Levy flight (7 < 1) to Levy walk 
(1 < 7 < 2). To our knowledge, this transition has been 
considered only as a notional possibility [l2j ]. 

In conclusion, we studied the stationary hole distri- 
bution p(x) produced by photo-excitation at an edge of 
n-InP slab and observed by broadsided luminescence at 
temperatures ranging from 300 K to 78 K. We discov- 
ered a giant increase of the hole spread in the sample 
over distances x exceeding 1 cm from the photo-excited 
edge. A power-law decline of the luminescence intensity 



4 



characteristic of Levy flight kinetics is observed with no 
change in the spectral shape. Our experiment provides a 
realization of the Levy flight of holes, in which the Levy 
distribution index 7 is controlled by the temperature and 
can be continuously lowered from 7 m 0.7 to 7 < 0.3 in a 
regular and well-understood manner. Such a control has 
never been available for formerly studied Levy processes. 
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